Superparabolic Level Glancing Models for Two-State Quantum Systems 
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Level crossing models for two-state quantum systems are applicable to a wide variety of physical 
problems. We address the special case of level glancing, i.e., when energy levels reach a degeneracy 
at a specific point of time, but never actually cross. The simplest model with such behaviour is 
the parabolic model, and its generalizations, which we call superparabolic models. We discuss their 
basic characteristics, complementing the previous work on the related nonlinear crossing models 
[Phys. Rev. A 59, 4580 (1999)]. 



I. INTRODUCTION 

The level crossing models provide an important tool for 
the description of nonadiabatic transitions. These tran- 
sitions occur when the energy levels of coupled quantum 
states are brought close to each other by driving them 
with external fields. The most prominent example of 
these models, the Landau-Zener (LZ) model, was intro- 
duced already in 1932 |lj-j3| in connection of atomic colli- 
sions (and by Majorana for atoms in time-varying mag- 
netic fields [4[), but because time-dependent level cross- 
ing problems have been studied over the years in con- 
nection with a wide variety of phenomena in physics and 
chemistry, for example with laser-atom interactions [5|- 
0] , laser-induced molecular dynamics H , slow l9Tfllt and 
cold [l2[ atomic collisions, molecular collisions [13j. neu- 
trino oscillations [HI, IH| and in Bose-Einstein condensa- 
tion [l6| , including evaporative cooling of atoms [ij} , out- 
couplers for atom lasers 1181- 211 ] and association of cold 
atoms into molecules |22l. |23|. In the recent years they 
have raised a lot of interest especially in nanophysics and 
quantum information processing as a way to coherently 
control qubits [lij and demonstrations of this have been 
done for example with solid-state artificial atoms (25l.[26|. 
Lastly, Landau-Zener type models have played an impor- 
tant role as a tool to understand the Kibble-Zurek theory 
of topological defect production and the dynamics of the 
quantum phase transitions [27-30]. 

First wc define some terminology for a quantum system 
with two coupled states and explicit time dependence in 
the Hamiltonian. The actual level crossing happens in 
the diabatic basis which consists of the bare states, i.e., 
system eigenstates when no coupling is present. The en- 
ergy difference between the diabatic states is termed as 
detuning, which becomes zero at the degeneracy point. 
However, such crossing appears only as an avoided cross- 
ing in the adiabatic basis (the basis of the instantaneous 
eigenstates of the time-dependent Hamiltonian) due to 
the coupling between the diabatic states. If the levels 
never reach degeneracy in the diabatic basis but get close 
doing so at some instant of time, transitions can happen 
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by tunnelling. In this paper we concentrate on a third 
kind of situation, which occurs as a limiting case between 
a proper crossing and a tunneling case, namely that the 
diabatic levels merely touch each other momentarily but 
do not actually cross. We call this a level glancing. 

In the original Landau-Zener model the diabatic en- 
ergy levels change linearly in time and the diabatic cou- 
pling is constant. Although this is a very crude assump- 
tion it has been applied very successfully over the years. 
The reason for this is that the nonadiabatic transition 
is located in the vicinity of the crossing point and one 
can usually linearize the diabatic levels in its neighbour- 
hood while the coupling does not vary much during this 
interval. In the recent years, however, there has been 
a growing interest to study more general dynamics than 
the one given by the single-crossing LZ case, for exam- 
ple in different interferometric schemes [3ll . |32| . One 
should note that in the original atomic collisions prob- 
lem with semiclassical trajectories, and which motivated 
Landau [l|], Zener |2| and Stiickelberg Q in 1932, the 
level crossing was traversed twice, and the phase differ- 
ence accumulated between crossings gives oscillations to 
the transition probabilities named after Stiickelberg, and 
they have been discussed and observed well before the re- 
cent interest in nanosystems, see e.g. 33]. In this article 
we follow the practice of refering to the single crossing 
case as the LZ model, and the "double-crossing" model 
with phase-related oscillations in transition probabilities 
as the LZS model. In collision physics the limit where 
the two crossing points in the LZS model approach each 
other relates to the case where the classical turning point 
overlaps with the region of degeneracy, and is accordingly 
considered as a breakdown of the LZS description. How- 
ever, that limit is also an example of a level glancing 
situation. Sadly, for historical reasons, the already more 
or less standard naming convention does not acknowledge 
Majorana's contribution Q. 

Another example of more general models are those 
with cubic-like detuning, i.e. the detuning is proportional 
an odd power of time. Such models have been analysed 
previously in Ref . [3"H , and they have recently raised some 
interest in the study of quantum phase transitions and 
adiabatic quantum computation 28]. In such cases one 
typically has to optimize between the computational time 
and the density of defects produced when crossing the 
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critical point, so dynamics different from the one given 
by the LZ model is needed, since in the adiabatic limit 
the transition time for the original LZ model increases 
exponentially with the coupling strength [35| . The cubic- 
like models deviate both quantitatively and qualitatively 
from the LZ model and were therefore dubbed essentially 
nonlinear crossing models in Ref. [34j . 

In this paper we consider the level glancing models that 
are reminiscent of the LZ and cubic-like models, having 
nonlinear parabolic-like time-dependent detunings and 
study the effect of these different time dependencies on 
the transition probability. One motivation behind this 
current work is the previously studied time-dependent 
parabolic model 136-38]. It can be used to describe sit- 
uations outside the scope of the usual LZ model, and 
it encompasses within a single model the cases of tun- 
nelling, double crossings and level glancing dynamics. 
Such parabolic time dependence of the diabatic energies 
has been recently applied e.g . in studies of laser-induced 
molecular dynamics [39L |40| . The parabolic model has 
a peculiar property: in the level glancing case the max- 
imum probability of a nonadiabatic transition is only a 
little bit over one half (3(|. Unlike in the LZ and cubic- 
like models, it does not reach unity in the adiabatic limit 
(actually, it goes to zero for the rather obvious reasons 
discussed later). Another important aspect is that the 
transition probabilities show oscillatory behavior for the 
level glancing situation as a function of relevant parame- 
ters, which is expected for any model with double cross- 
ings, but which are absent for the LZ model, for instance. 

We examine how the energy level dynamics affect the 
above properties by considering a set of level glancing 
models with superparabolic time dependence, i.e. the de- 
tuning depends on some even power of time. The os- 
cillatory character of the parameter dependence of the 
final transition probability can be understood on the ba- 
sis of the structure of the complex zeros of the adia- 
batic eigenenergies as explained by the Dykhne-Davis- 
Pechukas (DDP) theory (tlMi^. The results of the DDP 
theory are asymptotically exact in the adiabatic limit and 
although the structure of the zero points are rather simi- 
lar in both the cubic-like and superparabolic models, the 
diabatic limits of these models are completely different 
as mentioned before. 

The outline of the paper is as follows. In Sec. [TT] 
the basic equations and definitions are given. We also 
discuss the time evolution of the transition probabili- 
ties as well as complementary analytical approximation 
methods for the final transition probabilities, namely 
the above-mentioned Dykhne-Davis-Pechukas theory and 
various perturbation methods. In Sec. IIIII we present 
and analyze the results that were obtained by numeri- 
cal calculations and compare these to the approximative 
expressions. Finally, the discussion in Sec. IIVI ends the 
presentation. 



II. FORMALISM 

A. Basic equations 

We study the time evolution given by the Schrodinger 
equation (7j = 1) 

ij t m = H{t)m, (i) 

where H(t) is the general Hamiltonian of the two-level 
system with real- valued detuning A(t) and coupling 17(t), 
given explicitly by 

_ (A(t) n(t) \ . . 

H{t > - [n(t) -A(t) ) (2) 

and ij){t) = [ci(i),c 2 (i)] T , where ci(t) and c%(t) are the 
probability amplitudes of the diabatic basis states ip\ and 
ip2, respectively. As already mentioned, the level crossing 
in the diabatic basis is converted to an avoided one in the 
adiabatic basis, assuming that the coupling f2(t) does not 
vanish at the crossing. The adiabatic basis is formed 
by the instantaneous eigenstates and X2(t). The 

eigenvalues ±£ (t) of the Hamiltonian in Eq. ([2]) form 
now the adiabatic levels, with 

£(t) = v/tt(t) 2 + A 2 (t), (3) 

while the nonadiabatic coupling between the adiabatic 
states is 

j(t) = (xi(t)\X2(t)) = -{X2(t)\xi(t)) 

2(A(t) 2 + n{t) 2 ) ' 1 ' 

where the overhead dot stands for time derivative and 
one can fix the sign by fixing the relative sign of the 
basis vectors. 

B. Models 

We consider models where the coupling 17 (t) is con- 
stant and the detuning A(i) is directly proportional to 
some even power of time, 

17(t) = const. = 17 , A(t) = /3 N+ H N , (5) 

where (3 and 17o are real numbers and are both chosen to 
be positive and N — 2, 4, 6, ... is an even integer number. 
Now, instead of a crossing, the diabatic levels only touch 
each other at the point t = so we have a level glancing. 
An example of the energy level structures and couplings 
is depicted in Fig. [TJ These are fairly similar to other 
values of N as well. As the N increases, the parabolic- 
shaped energy levels transform to more rectangular ones 
and the two peaks of the adiabatic coupling get sharper. 
The corresponding physical picture is of course that we 
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drive the system fast to the resonance by changing the 
frequency of the constant-amplitude driving field, keep it 
there some time, and then the system is brought out of 
the resonance in a symmetrical fashion. 

Energy 





Time 



Coupling 



Time 



FIG. 1. The time-dependence of the energy levels for the 
superparabolic model with N = 4 and /3 = Q — 1. The 
diabatic levels and the corresponding coupling are drawn with 
solid lines while adiabatic ones are drawn with dashed lines. 



The case with N = 2 is called the parabolic model 
and it has been introduced in the context of atomic colli- 
sions [1, [l(| E3| and later used to study especially coher- 
ence effects related to multiple crossings and in situations 
where the LZ linearization fails [36l438| . Usually the di- 
abatic level energies in the parabolic model are given by 
at 2 — b so that the level glancing model ([5]) is actually 
only a special case of this with = 0. The parabolic 
level glancing model is able to give dynamics quite dif- 
ferent from the LZ model and is experimentally feasible, 
but one can not obtain with any parameters a higher 
probability of nonadiabatic transition than just over one 
half. We refer to the models with AT = 4, 6, 8, ... as su- 
perparabolic and show that with these models the situa- 
tion is improved. The models excluded from here where 
N is an odd integer number have been previously stud- 
ied in Ref. 34|. Note also that the case with N = 1 
is the exactly solvable LZ model. No exact solutions in 



a sufficiently simple or closed form seem to exist for the 
parabolic or superparabolic models, nor for any model for 
which the dynamics of the energy levels are even approx- 
imately given by Eq. §5§ near t — 0. This is in contrast to 
e.g. the Demkov-Kunike models with A(t) oc tanh(/3t); 
they are exactly solvable and reduce to the LZ model 
with certain parameter limits @, 13 . The models can be 
solved exactly, though, using the Zhu-Nakamura method 
based on Stokes parameters [IJ S3, EE]; but the re- 



sults are quite complicated, involving series solutions that 
eventually have to be evaluated using approximations 
and numerics (see, e.g., discussion in Ref. [46]). Our mo- 
tivation is to find reasonably accurate solutions, which 
have sufficiently simple expressions, and to examine in 
detail the DDP approximation. 

We take the system to be initially in state 2 so that the 
initial conditions are C2(— oo) = 1 and ci(— oo) = and 
we are interested in the asymptotic transition probability 
P. In the diabatic basis this is given by P — |ci(oo)| 2 . 
However, as the diabatic and adiabatic bases differ only 
by a phase factor as t — > ±oo and, unlike in the odd-iV 
cases, the basis states do not swap their labels this is the 
same as the probability P of a nonadiabatic transition. 
It should be also noted that now both the adiabatic and 
diabatic/sudden approximations correspond to P sa 0. 



C. Time evolution of populations 

Before discussing the final transition probabilities and 
their parameter dependence, we study the characteris- 
tics of the real time evolution of the populations by 
transforming to the rotating frame and using the Bloch 
sphere representation. The change of basis to the ro- 
tating frame is obtained by the unitary transformation 
17(f) = exp[-i(f>(t)a z }. Defining ij>(t) = U(t)ip(i) the 
Schrodinger equation (JlJ in this new basis reads 



A(t) - fat)) 



m. (6) 



The above change of basis shows that the time- 
dependence of the detuning can arise from the coupling, 
which is the case e.g. when couplin g at omic or molecular 
states with chirped laser pulses [8J, |40j, |49j . Generally, if 
we choose the time-dependence for the transformation to 
be 4> (t) = J A (s) ds, we will end up with the Hamilto- 
nian 



H(t) 



Q (t) e - 2 ^(*) 

ft (t) e+ 2i *(*) 



(7) 



This shows that any detuning can be alternatively 
mapped into a property of the coupling. Of course, as 
suggested by the Hamiltonian in Eq. ([6]), there can be 
time-dependence in both the detuning and the phase of 
the coupling, as in the so called winding models [Hof . 

Now we study the dynamics of the density matrix p in 
this basis given by the Liouville equation 



ip(t) = [H(t),p(t)}. 



(8) 



Furthermore, we define the three-component Bloch vec- 
tor R(t) = [x(t),y(t),z(t)] T , with 



X = P21 + Pl2, 

y = i(pi2 - P21) 

Z = PU — P22- 



(9) 
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For closed systems pxx + pn = 1 so that the third com- 
ponent of the Bloch vector gives also the evolution of 
the transition probability |ci(t)| 2 = |[l + z(i)]. From 
Eq. (JSJ we get the set of differential equations 

x = 20 sin (2(f>) z, 

y = -20 cos (20) z, (10) 
z = 20 cos (20) y - 20 sin (20) x. 

Equation (flO|) is equivalent to R = B x R. where the 
external field vector is given by 

B{t) = 20 (i) [cos (20(i)) , sin (20(t)) , 0] T . (11) 

The magnitude of the field vector is twice the value of 
the diabatic coupling, \B\ = 20 (t), and it is confined in 
the x-y plane with its direction determined by the angle 
20(i). 

In the case of the superparabolic models §5§ studied in 
this paper, the field vector traces a circle of radius 2O 
in the x-y plane. One can see here the difference to the 
usual Bloch vector representation given in the diabatic 
basis. There the field vector associated with the model 



in Eq. © is given as B = [O , 0, (3 N+1 t N ] T , so initially 
(t = — oo), it is aligned with the z-axis and its length is 
infinite. It is also confined in the x-z plane at all times. 
If the system is initially prepared to an eigenstate of the 
Hamiltonian, that is z(— oo) = ±1, then the Bloch vector 
will at first follow adiabatically the field vector but will 
eventually start to deviate and to precess around it as the 
direction of B changes. In the rotating frame the same 
initial conditions mean that the Bloch vector and the 
field vector are initially orthogonal and the adiabaticity 
can be understood from the fact that when the system is 
far away from the resonance point t — 0, the field vector 
is rotating so rapidly that its effects are averaged out as 
can be seen the particular choice of 4>(t) in Eq. The 
difference between the detuning functions of the super- 
parabolic level glan cing models and the cubic-like level 
crossing models in [34J is reflected to the behaviour of 
the field vectors, so that, in the rotating frame, the field 
vector will reverse its direction of rotation at t = for 
superparabolic models whereas for supercubic models it 
does not. 

By taking the initial conditions to be z(— oo) = — 1, 
y(—oo) = x(—oo) = and integrating the equations of 
motion (flUl) for the Bloch vector components, we get an 
expression for the time evolution of the z-component as 



z(t) = -1-4 I dt 0(t ){cos[20(i )] / dt ft(t )cos[20(£ )]z(t )+sin[20(i)] / dt 0(i )sin[20(i )}z(t )} 

J — OC 

dt"n(t)D,(t") cos[20(t') - 20(t")]z(t"). 
But we can also consider the one step before Eq. (fT2"|) . namely 



-1-4 



dt 

— oo ^ — OC 



(12) 



z(t) = -4O(i){cos[20(*)] / dt ft(t ) cos 20(t ) z(t ) + sin[2</>(£)] / dt 0(i ) sin[20(t )}z(t )}. (13) 



Multiplying both sides of Eq. (TT2)) by z(t) and defining 
the functions 

F c (t)= I dt'n(t')cos[2<j)(t)]z(t), 

T ( 14 ) 



Fs(t) 



dt 0(i )sin[20(t )]z(t ), 



we get the form 

z(t)z(t) = -4{F c (t)F c (t) + F s (t)F s (t)}. (15) 

By integrating Eq. (|15[) from — oo to t and imposing the 
initial condition z 2 (— oo) = 1 we get 



z'(t) = l-4{F?(t)+F s \t)} 



Let us now explicitly consider the superparabolic models. 
To this end, it is useful to first scale the time as r = (3t, 
leaving all the parameter dependence of the model Hamil- 
tonian given in the diabatic representation by Eq. ([5]) on 
the off-diagonal elements. We further define this new 
diabatic coupling as 



a = O //3 



(17) 



so that the limit a —> oo is now the adiabatic limit 
and a — > is the diabatic limit. By choosing the time- 
dependence of the transformation to the rotating frame 
(16) as 0(t) = t n+1 /(N + 1), we get Eq. ^ into the form 



5 



,2/(JV+l) 



z 2 (r) = I -a 2 [2 N (N+1)] 

In the limit of weak coupling, this simplifies to 
z 2 (r) = l-a 2 [2 N (N+l)] 2/iN+1 



dscos(s N+1 )z(s) ) + 



dssm(s N+1 )z{s) 



(C n+1 (t) + C N+1 (oo)) 2 + (S n+ i(t) - S N+1 (oo)) 2 



(18) 



(19) 



where 



C„(t) = / dscos(s n ), 
Jo 

S n (r) = f d SS m(s n ), 
Jo 



(20) 



r 



and n > 1. The limits of the integrals are given by 



r n /n) 

CJoo) = 1 1 ' cos (7r/2n) , 
n 

r(l/n) . 
o n (oo) = sin (Tr/Zn) , 



are the generalized Fresnel integrals [5l|, defined for r > so, in the weak coupling case, we get 
I 



(21) 



z 2 (oo) = 1 - 4a 2 



N + l 



2JV/(JV+1) 



r 2 



i 



N + l 



2{N+l) 



(22) 



Furthermore, as the assumption of weak coupling means 
that the diabatic coupling is too weak to change the pop- 
ulations considerably, it follows from the initial condition 
z(— oo) = —1 that z(r) ~ —1 for all time, so we can ap- 
proximate z 2 (oo) — 1 w — 2(z(oo) + 1). Now, from the 

I 



I 

definition of the Bloch vector in Eq. (j9]) we also have 
P = (1 + z(oo))/2, so that we finally get an expression 
for the final transition probability in the weak coupling 
approximation as 



P, 



pert 



N+l 



2N/(N+1) 



r 2 



i 



N + l 



2(N+1) 



(23) 



This approximation is expected to give accurate results 
only for very small values of the diabatic coupling and, 
indeed, it gives values over unity already with parameter 
values of a rs 1/2. To expand the range of validity of the 
perturbative approximation and to take the oscillatory 
character of the final transition probability into account 
one can also consider the Magnus approximation [52l l53l ] 
which can be obtained directly from Eq. (|23l) with the 
formula 



PMagnus — s hl 2 (\/ PBornj 



(24) 



which, of course, is always less than or equal to unity. Al- 
though this is an improvement for the models with larger 
N , in general this approximation is only of limited value. 
These approximations, along with the DDP approxima- 
tion which is derived in the next section, are compared 
to numerical results. 



D. Dykhne-Davis-Pechukas formula 



The Dykhne-Davis-Pechukas (DDP) formula is given 



by 



where 



p _ e -2ImD(t c 



D(t) = 2 [ £(s)ds, 
Jo 

and t c is defined by the equation 

£(t c ) = 0. 



(25) 



(26) 



(27) 



Equation (1231) gives the probability of nonadiabatic tran- 
sitions that is asymptotically exact in the adiabatic limit. 
This was proven by Davis and Pechukas [42j, who fol- 
lowed the original idea of Dykhne [4l| . The main assump- 
tions behind this formula are that £ (t) is non- vanishing 
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N = 2 



E 




FIG. 2. The complex zero points of the quasienergies are 
shown as black dots and the structure for the Stokes lines as 
solid lines for superparabolic models N = 2 and N = 4. We 
have chosen a — 1 so that the zero points lie on unit circle. 
The corresponding figures for higher values of N are similar, 
and in general there are three Stokes lines emanating from 
each zero point, and the points with equal imaginary part are 
linked by a Stokes line. 



ized to include the contributions of all the complex zeros 
relevant to the problem, as suggested in Refs. [42l. l53l. [56jj . 
so that 



P 



N 



k=l 



where 



T fe = 4z lim (t-t*) 7 (i), 



(29) 



(30) 



and j(t) is the nonadiabatic coupling defined in Eq. (Q}. 
From this formula it can be seen that the existence of 
multiple zero points lead to oscillations in the final state 
populations as the parameters are varied. 

The obvious problem in this case is to decide which of 
the points should be taken into account. As was shown 
rigorously by Joye and coworkers in Refs. [HI, [EH, this 
question is related to the global structure of the set of 
Stokes lines and the correct way is to include only the 
points connected by the Stokes line nearest to the real 
axis. However, they used an assumption that a limit- 
ing Hamiltonian exists as t — > ±oo which is not valid 
for many interesting models, including the present su- 
perparabolic models. Instead, it may be beneficial to 
include all the zero points in the Im(t) > plane in 
Eq. (p29"|) as studied in |34], [EH and this is also the pos- 
sibility we consider here. This viewpoint is supported 
by the fact that for the Demkov-Kunike models the sum 
over an infinite set of zero points actually produces the 
known exact result for all parameters [53j. 



1. Application to superparabolic models 



for real t (including t — ±oo), that t c is well separated 
from other zero points or possible singularities and that 
the Hamiltonian is analytic and single- valued at least in 
a region of complex t plane bounded by the real axis and 
the Stokes lines nearest to the real axis. The Stokes lines 
are defined as the level lines of D(t) with D(t) = D(t c ). 
Assuming that the zero points of £ (t) 2 are simple, it is 
easy to see by local analysis that there are three lines em- 
anating from t c with equal angles to each other as shown 
in Fig. [5] The zero points come in conjugate pairs, but we 
can restrict our considerations to upper half plane where 
Im(t) > 0. 

The DDP formula can generally be expected to a give 
a good approximation for P only near the adiabatic limit 
but for example for LZ model it gives the exact result, 
so it is correct for all the parameters. In the LZ case the 
only zero point on the upper half of the complex plane 
lies on the imaginary axis and Eq. (|2"5|) gives 

P = e- 7ra2 . (28) 

In the case that there exists more than one complex zero 
point, as is the case for the superparabolic models defined 
by Eq. ([5]) , the DDP formula (|25j) can be directly general- 



The zero points of the eigenvalues of the superparabolic 
Hamiltonian are 

T k = Q i/iv e w(2fe-i)/2iv 5 k = 1:2 ,...,N, (31) 

so the zero points lie on a circle of radius a 1 ^ . An 
example of the zero points and the structure of the Stokes 
lines is given in the Fig. [2j We get from Eq. (|26|) 

D(T k c ) = 7?e -(2fe-l)/2iV ! (32) 

where 

V = 2*W JV + 1 )/ JV (33) 

and 

where B(x,y) is the beta function [57j . The function i>n 
tends to unity as N increases. The factors in Eq. (I3TJ1) arc 
given by Tk = (— l) fe and the points can be grouped 
into pairs with the same imaginary part and opposite 
real part, so that the generalized DDP formula for the 
superparabolic models can be written in the form 



7 



DDP 



Ef/*(-l)V 



rjsin w 



(2fc-l)] 



^cos— (2fc - 1; 



(35) 



In many ways this is similar as the corresponding formula for the model in Eq. ([5]) with odd TV obtained in Ref. [34| , 



Tjodd _ » 
r DDP — 4 



UN-l)/2,_ 1)k VS in[^(2k-l)] 
/c— 1 V / 



^os^(2fc-l)]+(-l)<" +1 >/^ 



(36) 



except that with even TV the purely imaginary zero point 
is missing, which gives the last term in Eq. (1361) . Such a 
purely imaginary zero point is present in the LZ model as 
the only zero point, but in Eq. (|3"6")1 the largest contribu- 
tion in the large-a limit comes in fact from the first term 
of the sum and the role of the LZ-like term is suppressed. 
Thus the values for P show oscillations as a function of a 
also for odd values of TV, as long as TV > 1. This is a sim- 
ple mathematical explanation why the transition prob- 
abilities show oscillations that one would expect mainly 
for models that have interferometric character, i.e., two 
or more actual crossing or glancing points in the dia- 
batic basis. There is no obvious physical reason for such 
behaviour, although for even values of TV one can un- 
derstand the level glancing situation as a limiting case 
for the more general superparabolic model with double 
crossings [36l.l53j. 



III. COMPARISON OF METHODS 

We have applied the perturbative methods, the DDP 
methods and numerical evaluation to the superparabolic 
level glancing models. The asymptotic probabilities 
Pd d p and PMagnus for the parabolic and superparabolic 
models are shown in Figs. [3] and 0] along with the results 
obtained from numerical calculations as the coupling a 
is varied from to 6. They are oscillating functions with 
respect to a in contrast to the monotonous relationship 
between the probability of transition in the LZ model 
and its adiabatic parameter. In Fig. @]we show the same 
curves but with logarithmic scale for the coupling in order 
to study the approximations near the maximum values 
of the transition probability. We demonstrate the diver- 
gence of the weak coupling approximation Pp er t in this 
figure. But before making the full comparison, let us 
discuss a few special points in light of the shown results. 

One can see that the validity of the generalized DDP 
approximation extends to increasingly smaller values of a 
as TV grows and apart from small phase shift it is in a good 
agreement with the numerical result for all values of a. 
Let us next consider the behaviour of the superparabolic 
models as TV — » oo. As already mentioned, the behaviour 
of the diabatic energy levels in the limit of large TV indi- 
cates that one could expect the time-evolution of the sys- 
tem to be divided into two parts, i.e., to a non-resonant 
part where the time-evolution is not affected by the cou- 
pling, and to a resonant period where the dynamics is de- 
scribed by the well-known Rabi model. The Rabi model 



corresponding to the superparabolic model in this reso- 
nant period consists of a constant rectangular pulse, of 
magnitude a, which couples the levels. Thus one expects 
to obtain the final transition probability P by just deter- 
mining the duration of the resonant period and inserting 
it to the well-know expression for resonant Rabi oscilla- 
tions. To get an estimate for this duration time we use 
the same argument that is generally used e.g. in esti- 
mating the dynamically relevant time region for the LZ 
model (this crude approximation is sufficient if a is not 
too large, see Ref. [35|]). Then the start and the end of 
the transition region is determined simply from the con- 
dition A(t) = O(r), which gives At = 2a 1 / JV , so that in 
the limit TV — s- oo the final transition probability accord- 
ing to the Rabi formula is 



Ppabi = sin 2 [2a] 



(37) 



This is consistent with the results of the DDP approx- 
imation with large TV as can be seen from Fig. \5\ As 
expected, deviations begin to occur when a increases as 
then our crude estimate for resonant period begins to 
falter. 

To evaluate the difference between the generalized 
DDP approach and the original DDP theory, we con- 
sider Pddp when one includes in the sum only the zero 
points closest to the real axis. We denote this original 
theory prediction by P^bp- As it has been shown there 
are two such points located symmetrically with respect 
to imaginary axis and those correspond to the first term 
in the summation in Eq. (|35[) . so that 



P«^ = 4e- 2 " sin [w] sin 2 



r\ cos 



2TV 



(38) 



In the case of weak coupling, a ~ 0, this becomes 



P. 



(i) 



DDP 



(2N+1)/N ^[(2N+l)/(2N)} (JL) 
a r 2 [(3TV + l)/(2TV)] \2NJ 



which is of interesting form when compared to the proper 
weak coupling limit in Eq. ([23| . When TV increases they 
both tend to limit a — > as a 2 , but for finite TV the 
two expressions are slightly different. For the generalized 
DDP result in Eq. (|35l> . we get a sum with all terms 
proportional to a 2 in the small a limit, but the sum 
does not appear to have any simple closed form for easy 
comparison. 

Let us now proceed with the general examinations of 
the results. We can see from the Fig. [5] and Table Q] 
that the maximum probability is, as we expect, enhanced 
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FIG. 3. The probability P of a nonadiabatic transition is 
plotted as a function of the coupling a for the values N — 2, 6 
and 10. The dots indicate the numerical solution and the solid 
line is the generalized DDP result. 



N 


p 

J max 


Q 


2 


0.577 


0.68 


4 


0.792 


0.68 


6 


0.885 


0.68 


8 


0.929 


0.69 


10 


0.952 


0.70 


12 


0.965 


0.70 


14 


0.974 


0.71 


16 


0.981 


0.71 



TABLE I. The maximum values of P that appear in Fig. [6] 
and also the values of coupling cto with which they can be 
obtained. 

by increasing N, though the rate of the increment de- 
creases steadily. At the same time, the value of coupling 
cto needed to obtain the maximum value increases only 
slowly. 

It seems that the parameter value cto f° r which the 
maximum transition probability is obtained is somewhere 
in between the region of validity of the perturbation and 
DDP approximations. It is, however, evident from Figs. [3] 



FIG. 4. This figure is similar to Fig. © but the vertical scale 
is normal whereas the horizontal scale is logarithmic. In this 
figure, the thick solid line is the numerical result and the thin 
line is the DDP result. Furthermore, the dot-dashed line in 
this figure represents the weak coupling approximation Pp er t 
and the dashed line is the Magnus approximation. 



and |4] that the generalized DDP method is in good agree- 
ment with the numerical results for the superparabolic 
models in a large part of the parameter region and not 
just in the adiabatic limit, and its accuracy gets better as 
N increases. Because of the oscillating character of the 
transition probability, small error in the phases can lead 
to notable deviation from the true value, and there in- 
deed is some phase difference in the DDP and numerical 
results but the generalized DDP results still catches all 
the essentials of the numerical results. For small values 
of a, the perturbation approximation is more accurate, 
its region of validity remaining somewhat constant. The 
Magnus approximation improves the perturbation result 
further but in a limited fashion. 

It is clear that to get a quantitatively accurate approx- 
imation using the generalized DDP method, one has to 
include all complex zero points into the expression (|29|) . 
Taking into account only the pair of complex zero points 
closest to the real axis and which are connected by the 
closest Stokes line gives qualitatively somewhat correct 
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N= 80 




FIG. 5. This figure demonstrates the behaviour of Pddp 
(solid line) and Pnabi (dashed) for large N. Here we have 
chosen N = 80. 



1.0 
0.9 
0.8 
0.7 
0.6 



2 4 6 



10 12 14 16 



N 



FIG. 6. The maximum probabilities that can be obtained for 
each model as a function of N. 



behaviour and for example the parameter value for ob- 
taining the maximum probability but also values over 
unity. 



IV. CONCLUSIONS 

Although the Landau-Zener model has become, for 
many reasons, a paradigm in studies of coupled time de- 
pendent quantum states, the parabolic model and espe- 
cially the level glancing aspect have not been considered 
in detail, and the superparabolic models have not been 
considered at all, to the best of our knowledge. A level 
glancing situation is, of course, much less likely to occur 
in nature, although increasing interest and development 
of tools for control and engineering of quantum states 
will make it likely that such models can be tested and 
applied in experimental physics. Our motivation for the 
present study has been three-folded. 

Firstly, we wanted to understand the nature of the level 
glancing dynamics and especially the oscillating charac- 
ter of the transition probabilities, i.e., how the proper- 
ties of the parabolic model extend to the superparabolic 
models. We have shown that the superparabolic models 
provide dynamics that follows qualitatively the parabolic 



model, and eventually starts to resemble the Rabi model, 
with a sudden approach to the resonance, steady resonant 
Rabi oscillations, and then the freezing of the value of the 
transition probability to the moment of the sudden move 
away from the resonance. 

Secondly, we have complemented the picture of essen- 
tially nonlinear models as studied in Ref. 34], by intro- 
ducing the superparabolic models, and showing how the 
structure of the complex zeros reflects to the transition 
probability. It is clear that the Landau-Zener model is in 
many ways unique in its monotonous change of transition 
probability as a function of coupling strength, whereas 
both the cubic-like and the superparabolic models will 
display oscillations, and the superparabolic models allow 
one to reach full transition only asymptotically. 

Thirdly, the improvement of the adiabatic solution 
(DDP) by addition of all zeros of the eigenenergies in 
the complex plane lacks rigorous proof so far, and its 
practicality must be evaluated for each model indepen- 
dently 53] . Our work shows that for the superparabolic 
level glancing models the inclusion of all such zeros in 
the upper half of the complex plane improves the qual- 
ity of the solution and extends its usefulness beyond the 
adiabatic limit. 

One aspect of the Landau-Zener model is the intrin- 
sic phase difference between the quantum state ampli- 
tudes which the level crossing transition dynamics pro- 
vides Q. This becomes visible in transition probabil- 
ities for any double or multiple crossings situation, al- 
though the dynamical phase evolution of the amplitudes 
between the crossings tends to dominate any oscilla- 
tions (36( | . The other models, including parabolic and su- 
perparabolic models, give rise to a similar phase factor. 
The Landau-Zener model has been applied to periodi- 
cally driven systems [58l-l60| , and it would be interesting 
to see if a system with periodically occurring level glanc- 
ing events gives similar dynamics. A second aspect to 
consider is the role of noise on both the cubic-like as well 
as parabolic and superparabolic models. The Landau- 
Zener model as such is not very much affected by the 
noise [6l|, but any double or multiple crossing situation 
will be strongly affected by phase degradation [25| , offer- 
ing, on the other hand, a tool for analysing decoherence of 
quantum states. Control of two-state systems as well as 
their decoherence is important for quantum information 
and quantum computing [62j . The essentially nonlinear 
models all show oscillations in the transition probabili- 
ties, even though there is only a single level crossing or a 
level glancing event. Their sensitivity to decoherence is 
an interesting and open question, and one of the topics 
for further studies on level glancing models. 
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